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Quantum search using non-Hermitian adiabatic evolution 
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We propose a non-Hermitian quantum annealing algorithm which can be useful for solving complex 
optimization problems. We demonstrate our approach on Grover's problem of finding a marked item 
inside of unsorted database. We show that the energy gap between the ground and excited states 
depends on the relaxation parameters, and is not exponentially small. This allows a significant 
reduction of the searching time. We discuss the relations between the probabilities of finding the 
ground state and the survival of a quantum computer in a dissipative environment. 
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Many physical and combinatorial problems 
associated with complex networks of interact- 
ing degrees of freedom can be mapped to equiv- 
alent problems of finding the minimum of cost 
function or the ground state of a corresponding 
quantum Hamiltonian, Ho, [THi)]- One of the 
approaches to find the ground state of Ho is 
quantum annealing (QA) which can be formu- 
lated as follows. Consider the time-dependent 
Hamiltonian H(t) = Ho + T(t)Hi, where Ho is 
the Hamiltonian to be optimized, Hi is an aux- 
iliary "initial" Hamiltonian, and [Ho, Hi] 7^ 0. 
The coefficient, T(i), is a control parameter, 
and r(t) decreases from very high value to zero 
during the evolution. 

One starts with the ground state of Hi as 
the initial state, and if T(t) is slowly decreasing, 
the adiabatic theorem guarantees approaching 
the ground state of Ho, at the end of the com- 
putation, assuming that there are no energy 
level crossings between the ground and excited 
states. So, the quantum optimization algo- 
rithms require the presence of a gap between 
the ground state and first excited state. How- 
ever, in typical cases the minimal gap, g m , 
is exponentially small. For instance, in the 
commonly used quantum optimization n-qubit 
models, the estimation of the minimal energy 
gap yields: g m « 2""/ 2 0, |, EMI- This 
increases drastically the total computational 
time, and from a practical point of view the 
advantage of the method is lost. 

Recently (i~3j . we have proposed a non- 
Hermitian adiabatic quantum optimization 
with the non-Hermitian auxiliary Hamiltonian. 



We have shown that the non-Hermitian quan- 
tum annealing (NQA) provides an effective 
level repulsion for the total Hamiltonian. This 
effect enables us to develop an adiabatic the- 
ory without the usual gap condition and to de- 
termine the low lying states of Ho, including 
the ground state. Some interesting suggestions 
for implementation of non-Hermitian architec- 
tures by realization of "Ising machine" based 
on mutually injection-locked laser systems were 
recently discussed in [Til . [H^ . 

In this Letter, we apply the NQA to Grover's 
problem (l6| . i.e. finding a marked item in an 
unstructured database. 

Consider a set of N — 2™ unsorted items 
among which one item is marked. The related 
Hilbert space is of dimension N. In this space, 
the basis states are written as \i) (i=l,2,. . . ,N), 
and the marked state is denotes as \m). The 
task is to find the marked item as rapidly as 
possible. 

The Hamiltonian whose ground state is to 
be found, can be written as: Ho = — \m)(m\. 
Its ground state, marked as |m), is unknown. 
The auxiliary Hamiltonian is given by "Hi = 

-IVoXM where |Vo) = (l/V^)E£i V) is its 
ground state with energy Ef = — 1. For both 
Hamiltonians, Ho and Hi, the rest of eigen- 
states have the N — 1-times degenerate energy 
E r = (r = 2, 3, . . . , N). (Our choice of the 
Hamiltonian is different from the Hamiltonian 
considered in refs. (l7l - [20j by a total shift on 
the unit matrix.) The total time-dependent 
non-Hermitian Hamiltonian is chosen as fol- 
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lows: H T (t) =% + h{t)%x, where 

«*>-{?-,& "* ,ir o 

We denote 7 = (5 + i<5) /r, where 5 and 5 are 
real. In what follows we assume that 8 <C 1. 

The adiabatic quantum search algorithm 
consists of (i) preparing the system in the initial 
state, \ip(0)) = \ipo), and (ii) performing an evo- 
lution by applying the non-Hermitian Hamilto- 
nian, 7i T (t), during a time, r. At the end of 
evolution, the non-Hermitian part of the total 
Hamiltonian disappears. Then, if the evolution 
is sufficiently slow, the system is remained in 
its ground state, which will be the ground state 
of the Hermitian Hamiltonian, Hq. 

We start with the solution of the eigenvalue 
problem for H T (t). This yields (N-2)-times de- 
generate highest eigenvalue, E2 = 0, and two 
lowest eigenvalues, Eq and E\, which are given 

by 



Choosing the orthonormal basis as 



(2) 
(3) 



where Q(t) = 
e(t) = h(t) + 1 



2h(t) cos a + 1 and 
We set sin(a/2) = 1/y/N. 
The energy gap between the ground state and 
the first excited state is given by |A25(t)| = 
\y/h 2 (t)-2h(t) cosa + l|. For N > 1/6 one 
can show that the minimum of the energy gap 
is given by |A£| min = S/^fTS 2 + 0{1/N). 

In the two-dimensional subspace spanned 
by the vectors, \ipo) and |m), we choose 
an orthonormal basis as |V>o) and \ij)x) = 
(sin(a/2)|-0 o ) — |m))/cos(a/2). We comple- 
ment it to the basis of the iV-dimensional 
Hilbert space by adding (N — 2) vectors \ipk) 
(k = 2, . . . , N — 1), which form the orthonor- 
mal basis of the orthogonal (N — 2)-dimensional 
Hilbert subspace. Then, an arbitrary state, 
can be expanded as |*(*)) = c (t)\ip ) + 

Inserting this expansion into the Shrodinger 
equation, id/dt\V) = H T |*(*)>, we find that 
the differential equations for the coefficients, 
Co(t) and ci(i), do not involve the coefficients, 
c k (t) (k = 2, . . . ,N - 1). Then, effectively the 
iV-dimensional problem is exactly reduced to 
the two-dimensional one. So, it is suffices to 
confine our attention to the two-dimensional 
subspace. 



• 0/ 

l j,\ipi) = (o)}> one can wr it e the cor- 
responding effective (non-Hermitian) Hamilto- 
nian as 



sit) n(t) 

rief{t) = Y ~2~ ' 



(4) 



where fi(t) = (sin a, 0, h(i)— cos a) is a complex 
vector, and <x denotes the Pauli matrices. 

We denote the (right) instantaneous eigen- 
vectors, corresponding to the eigenvalues, 
E a (t), as \u a (t)) (a — 0,1). One can show 
that K(0)) = |Vo) +0{1/N), and \u (t)) -> 
+0(1/N), as t — )■ t. 

For the two-level system (TLS) governed by 
the effective non-Hermitian Hamiltonian (|4]), 
the wave- function can be written as \tp(t)) = 
Co(i)IV'o) + ci(t)\ipi). We assume that the evo- 
lution of TLS starts at to — in the state 
\ip(0)) — \ipo). This implies the following ini- 
tial conditions: Co(0) = 1 and Ci(0) = 0. 

Writing c a = U a (t) exp (| J" ' e(t)dt) , and em- 
ploying the Schrodinger equation for the TLS 
governed by the effective Hamiltonian, T-L e f{t), 
we obtain the Weber equation for the new func- 
tions U a (z), 



d 2 , . , 1 z 2 



iv)u ,i(z) = 0, 



(5) 



where z{t) = e l7T / 4 (-y(T — t) — cosa)/^ and 
v = sin 2 a/47. 

Solutions of Weber's equation are given by 
the parabolic cylinder functions, D-i V {±z), 

U (z) = {AD_ iv {z) + BD^(-z)), (6) 

U x (z) = Vw(BB-„-i(-z) - AD-i V -i(z)). 

(7) 

The constants, A and B, being determined 
from the initial conditions, are found to be: 
A = D_ w _i(-z )F(l + iv)/y/5nw and B = 
D-i V -i(zo)F(l + iv)l\j2nv 1 where we set zo = 
z(0). " 

It is assumed that the quantum measurement 
will determine the state of the quantum sys- 
tem at t > t. We denote the final state of the 
system as \ip T )- Then, the probability, P n , of 
finding the system in a given state, \n), can be 
written as, 



Pn = 



IWr)| ; 



2 ' 



(8) 
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This yields the (intrinsic) probability of transi- 
tion ) — > \tpi) as 



Pr(t) = 



co(t)l 2 + |cx(t)j= 



(9) 



Thus, P t {t) is the probability of the system 
being in the ground state at the end of the evo- 
lution. 

Using the functions C/o,i(z), we write the 
probability of transition, P r (r), as 



PAr) 



1 



1 



WTWW 



(10) 



e 1 ™/ 4 cos a/^Fj, and for 
" /4 /V7- 



where z T — z(r) = 
N ^> 1 we obtain: z r w — e 

To estimate the probability of transition, we 
apply asymptotic formulas for the parabolic 
functions. This yields 



U (z T ) _ e -W2 e -4/2r(l + j v ) 



/27TZq 



(11) 



Inserting (fTT|) into E q. El) , we obtain the 
Landau-Zener formula [2ll . |22j for the Hermi- 
tian quantum search (5 = 0), 



Pr(r) = 1 



-2-kv 



(12) 



where v = (r/gN). We conclude that P t (t) ~ 
1) if T > 5 AT. Thus, to obtain the probabil- 
ity close to 1 to remain in the ground state at 
the end of evolution, the computational time 
should be of order N. In fact, this result is 
equivalent to the well-known result on the com- 
plexity of order N provided by quantum adia- 
batic evolution approach [18[ , which is the same 
as in the classical search algorithm. 

For the NQA, assuming N 3> 1, we ob- 
tain the following rough estimate of the com- 
putational time: r > (g 2 /5) In N. Thus, the 
non-Hermitian quantum search has complex- 
ity of order In N, which is much better than 
the quantum Hermitian (global) adiabatic algo- 
rithm. Also, this complexity is certainly better 
than one of the adiabatic local search algorithm 
that has total running time of order y/~N 17 1. 



In Fig. [T]we present the results of our numerical 
simulation. For the Hermitian QA (S — 0) the 
transition probability is: P T « 3 • 10~ 8 ; and for 
the NQA with weak dissipation, 5 — 0.0025, the 
transition probability is: P T = 1 (r = 1.5 • 10 4 ). 

Nonlinear NQA. - We define the survival 
probability of the lossy system as the trace of 
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FIG. 1: Left panel: The transition probability, 
P T as a function of the scaled time, s — t/r 
(8 — 0). Right panel: The survival probability, 
Ps{s) (dotted blue line), and the transition proba- 
bility, P T (s) (red line) (S = 0.0025). In all cases: 
q = 2,r = 1.5 ■ 10 4 ,iV = 2 40 . 
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FIG. 2: Nonlinear NQA. The transition probabil- 
ity (left panel), and the survival probability (right 
panel) as the functions of the scaled time, s = t/r 
(g = 2,5= HT 4 , t = 5 ■ 10 4 , N = 2 40 ). 



the density matrix, P s (t) = Trp(t). Using the 
asymptotic formulas for the Weber functions, 
one can show that for N 3> 1, the asymptotic 
behavior of the survival probability is given by: 
P s (t) w e~ st . (See Fig. Q]) Then, one can see 
that the conditions to obtain high probabilities 
for (i) finding the ground state, leading to in- 
equality, t > (g 2 /5) In N, and (ii) survival of 
qubits, St < 1, are not compatible. A compro- 
mise can be found by using a local adiabatic 
evolution approach [17]]. 

We rewrite the total time-dependent non- 
Hermitian Hamiltonian as, 



H T {t)=Ha + h (l-f(t))Hi 



(13) 



where ho = (g + i5), and fit) is a mono- 
tonic function of t. For concreteness, we choose 
g = 2, and impose the following boundary con- 
ditions: /(0) = and /(r) = 1, where r de- 
notes the computational time. 
We choose f(t) solution of, 



dt t \ 



(i - V) 



(14) 
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FIG. 3: Nonlinear NQA.The transition probabil- 
ity (left panel) and the survival probability (right 
panel) as functions of the scaled time, s = tjr 
(g = 2,S = 7.5 • KT 5 , r = 5.5 • 10 4 , N = 2 40 ). 



where (3 = arctan(l/<5). 
gration, we find 



Performing the inte- 



t= 1 + ^aretan i^il. (15) 
By inverting this function we obtain, 

/(t) = i + £tan(^(|-l)). (16) 

From here it follows that /(t) = 1, and the 
computation time is r. 

In Figs. [2] and [3] we present the results of nu- 
merical calculations for different choice of pa- 
rameters, S and r. Our results show that the 
nonlinear NQA can be realized with the tran- 
sition probabilities, P T « 1.2 ■ 10~ 2 and P s « 
1.6 • 10" 2 . The computation time, r rj 5.5 • 10 4 , 
is better than the time of quantum search pre- 
dicted by the Grover algorithm, r = y/N ~ 10 6 
(for n = 40). 

Conclusion.- The field of quantum adiabatic 
computation is well-established, and many use- 
ful results are discussed in the literature. One 
of the main problems of this approach is that 
the energy gap between the ground state to be 
found and the excited states is generally ex- 
ponentially small. This requires exponentially 
large computation times. On the other hand, 
in the dissipative (non-Hermitian) regime, the 
energy gap is defined by the relaxation parame- 
ters, and may not be exponentially small. (See 



also [lj,ll5(.) In this case, the computational 
time can be significantly reduced. But then, an- 
other problem appears-the quantum computer 
has a finite probability to be destroyed (which 
happens anyway). One way to overcome this 
problem was discussed in [14,15], where both 
dissipation and external pumping in the locked 
laser system was used to model the Ising sys- 
tem in its stationary ground state. But still 



many theoretical and experimental issues must 
be resolved in order to build this type of "Ising 
machine" . 

The results presented in our paper demon- 
strate that non-Hcrmitian quantum computa- 
tions can be used for two purposes. One is 
to use non-Hcrmitian quantum algorithms to- 
gether with the use of classical computer to re- 
duce computation time. We are in the process 
of demonstrating this option for some classes of 
Ising models. Another purpose is to build a real 
"non-Hermitian quantum computer" (NHQC) 
to solve specific problems rapidly. As was 
demonstrated in this paper, there will be a 
tradeoff between the probability of finding the 
desired outcome and the probability of survival 
of the computer. As our results show, there are 
useful ways to improve the performance of the 
NHQC. 
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